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In recent years there has been a growing interest in the study of shape formation using modern 
responsive materials that can be preprogrammed to undergo spatially inhomogeneous local defor¬ 
mations. In particular, nematic liquid crystalline solids offer exciting possibilities in this context. 
Considerable recent progress has been made in achieving a variety of shape transitions in thin sheets 
of nematic solids by engineering isolated points of concentrated Gaussian curvature using topolog¬ 
ical defects in the nematic director field across textured surfaces. In this paper, we consider ways 
of achieving shape transitions in thin sheets of nematic glass by generation of non-localised Gaus¬ 
sian curvature in the absence of topological defects in the director field. We show how one can 
blueprint any desired Gaussian curvature in a thin nematic sheet by controlling the nematic align¬ 
ment angle across the surface and highlight specific patterns which present feasible initial targets 
for experimental verification of the theory. 


I. INTRODUCTION 

Complex shape transitions driven by inhomogeneous 
local deformation patterns are ubiquitous in nature and 
are observed at a variety of length scales, ranging from 
the cell walls of plants and bacteria to macrostructures 
such as plant leaves [Tj . Examples of such local deforma¬ 
tions include differential growth in biological structures 
and nonuniform mechanical responses of biological tis¬ 
sues to external stimuli such as humidity [2]. Modern 
responsive materials can be preprogrammed to undergo 
spatially inhomogeneous expansions and contractions in 
response to external stimuli such as heat and light. Ex¬ 
amples include thermoresponsive hydrogels 33 and ne¬ 
matic liquid crystalline solids [4|. There has been con¬ 
siderable interest in recent years in the prospect of pre¬ 
programming desired shape transformations that can be 
remotely activated in such materials EHZ]. 

Liquid crystalline solids consist of long, flexible molec¬ 
ular chains that are sufficiently cross-linked to form a 
solid. At sufficiently high temperatures, the molecular 
directions are randomly distributed and the material is in 
the isotropic phase. When the temperature is below some 
critical value, the rod-like molecular elements become lo¬ 
cally aligned about the director n and the material is said 
to be in the nematic phase. Liquid crystalline solids expe¬ 
rience elongations and contractions in response to light, 
heat, pH, and other stimuli that change the molecular 
order. Of particular significance are nematic glasses 0 
and elastomers [9]. Both have spontaneous deformation 
tensors of the form 

F = (X— \~ v )n (gm + A - " Id 3 , (1) 

where M 3 denotes the identity operator on R 3 . That is, 
a local scaling by A along the director n and a scaling 
by A - " perpendicular to n is observed upon exposure to 
stimulus. The parameter v is known as the opto-thermal 
Poisson ratio and relates the perpendicular and parallel 
responses H;. 

In nematic glasses, molecular chain motion is highly 
limited by the cross-links and the director is not in¬ 


dependently mobile from the elastic matrix as it is in 
elastomers. For nematic glasses, we typically have A £ 
(0.96,1.04) and v £ (|, 2), while for elastomers A £ (|, 4) 
and v = A. In nematic glasses, the director field changes 
only through convection by mechanical deformations, 
which allows for a feasible patterning of the director field 
at the initial time of cross-linking and the subsequent 
guarantee that the pattern will not be modified by “soft 
elasticity” mediated by director rotation [9j. Consider¬ 
able control of director fields is achievable through a va¬ 
riety of techniques including the use of electric and mag¬ 
netic fields, surface anchoring, and holography. Thus, 
the matrial response in nematic glasses can be prepro¬ 
grammed by precisely setting up a desired director field 
pattern immediately before the glass is formed via cross- 
linking mi¬ 
lt has come to our attention that a recent ambitious pa¬ 
per by Aharoni et al. nm addresses the question of shape 
transformations that are achievable through the pattern¬ 
ing of general smooth director fields on flat sheets, as was 
the original motivation for much of this work. Nonethe¬ 
less, we hope that this paper will add to their contribu¬ 
tion by highlighting specific patterns which present fea¬ 
sible initial targets for experimental confirmations of the 
theory. We also introduce the notion of orthogonal du¬ 
ality of director fields and consider its implication for 
Gaussian curvature distributions. Nematic patterns on 
initially curved surfaces are considered towards the end 
of the paper. 

The value of A and v are typically dependent on the 
strength of the stimulus. In the following analysis, we 
assume that the shape transformation is achieved by a 
sudden exposure to stimulus of fixed strength, so that 
there is a sudden activation of the prescribed local de¬ 
formations for fixed values A and v. For definiteness, we 
assume that A < 1, so that the local deformations consist 
of a contraction along the director and an elongation in 
the perpendicular directions. We use the convention that 
Latin indices run over 1, 2, 3, whereas Greek indices take 
values 1 and 2 . 
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II. NON-EUCLIDEAN LOCAL 
DEFORMATIONS AND NEMATIC SURFACES 

A configuration of a body in R 3 can be specified by a 
smooth injective immersion $ : fl —> R 3 , where 0 C R 3 
is a domain. The configuration defines curvilinear coor¬ 
dinates x = (x l ) = (a; 1 , a; 2 , ir 3 ) G 0 for material points 
throughout the body. The map <1? : il —> R 3 induces a 
Riemannian metric g on O with covariant components 

g ij (x) = d i *-d j *. (2) 

The Christoffel symbols are defined in terms of the com¬ 
ponents of the metric tensor by 

I "\ji — ^{&jgu T digji digtj) (3) 

and T*- = g kl Tiji, where ( g = (g^) -1 . It is well-known 
that the components 


integrating over the shell thickness h. In the case of a 
homogeneous and isotropic elastic material, we obtain 
the reduced energy functional 

£ = / W(xi,x 2 )y/\d\dxidx 2 

J UJ 

= / (Ws(xi,i 2 ) + W B (ii,i 2 )) VW\ dxidx 2 (6) 

J oj 

where 

Ws(xi,x 2 ) = 1 ^A aM5 {a a p - a a0 )(a lS - d lS ) ( 7 ) 

Wb(i i,x 2 ) = (b a p - b a p)(b^ s - b lS ) (8) 

are the stretching and bending energy contributions, re¬ 
spectively P’2]- The reference fundamental forms d a p 
and b a p are related to the three-dimensional reference 
metric g via 


Rujk = djT iU — d k Tiji + rur fc j m — (4) 

of the Riemann curvature associated with any such in¬ 
duced metric uniformly vanish at all points in [Ill- 

Processes such as differential growth in biological 
structures can determine a reference geometry that is 
characterised by a reference Riemannian metric g that 
describes the prescribed rest distances associated with 
the underlying growth law. Typically, the underlying 
growth law gives rise to a reference metric that is non- 
Euclidean in the sense that its Riemann curvature tensor 
has non-vanishing components [5]. 

As the body deforms in response to the underlying lo¬ 
cal deformations, it must assume a configuration in three- 
dimensional Euclidean space. Thus, its realised configu¬ 
ration is characterised by a Riemannian metric g that 
is induced by a configuration $ : fl —»• R 3 . The met¬ 
ric g is referred to as the actual metric. If the under¬ 
lying local deformations are geometrically incompatible 
with Euclidean space, the prescribed rest distances de¬ 
scribed by g cannot be realised everywhere. This gener¬ 
ates a residual strain field throughout the body given by 
e(x) = \ ( g(x ) — g(x)). For a hyperelastic material, the 
elastic energy is given by a functional 

E[*\= f W(x,g(x))dV, (5) 

Jn 

where dV = \f\f \ dx\dx 2 dx 3 . At each point x, the energy 
density function W(x, g(x)) vanishes if and only if g(x) = 
g(x) P2]. 

Consider a standard domain of parametrisation of the 
form Q = ui x (—§,§) for a thin shell of thickness h , 
where ui corresponds to the mid-surface and the third 
coordinate £3 G (— f > f) measures the distance along the 
normal to the mid-surface. We can derive a reduced two- 
dimensional model formulated in terms of the first and 
second fundamental forms a a p, b a p of the mid-surface 
$ (cj) by using the Kirchhoff kinematic assumption and 


^ck/3 (*^1 1 *^ 2 ) 9ot(3 |ai3=0? b&fl (*^1 5 *^ 2 ) cy &39cx(3 1x3=0’ 

(9) 

The components of the elasticity tensor A a ^ lS [12] are 
given by 

A "' h ‘ = 4(1^2) M*** 

+ i(i-i'»L)(<s“ 7 «' , ' + a‘" i a»T)], (10) 

where the Young modulus Y and Poisson ratio v e \ are 
related to the Lame coefficients A e i and g e \ by 


Ve 1 

1 - t'el 


Aei 

A e i + 2 /r e i 


Y = 2/X e i (1 + I/ e i). 


( 11 ) 


As h —> 0, the equilibrium configuration manifests as 
an isometric immersion of the reference two-dimensional 
metric a a g. That is, 


a ct/3 = a a/3 as h —► 0. (12) 


More precisely, the analysis of the thin sheet limiting 
behaviour has been placed on a more rigorous mathe¬ 
matical foundation in recent years using the concept of 
T-convergence. In Friesecke et al. na the nonlinear 
bending theory of plates due to Kirchhoff is derived as 
the T-limit of the classical theory of 3D nonlinear elas¬ 
ticity, under the assumption that the classical 3D elas¬ 
tic energy per unit thickness h scales like h 2 . The non- 
Euclidean version of this result is derived in m under 
the same scaling law applied to the energy functional 
([5]), yielding a natural non-Euclidean generalization of 
the Kirchhoff model with a corresponding 2D bending 
energy functional. Necessary and sufficient conditions 
for the existence of a W 2,2 isometric immersion of a given 
2D metric into R 3 are also established. In particular, it is 
shown that if finite bending energy isometric immersions 
of the metric exist, then the minimizers of the 3D elastic 
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energy © converge in the vanishing thickness limit to an 
isometric immersion of the mid-surface metric, which is a 
global minimizer of the bending energy over all isometric 
immersions m- 

Consider a director field n : Q —> R 3 given by n = 
n(xi, X 2 , X 3 ) = rii(xi,X 2 ,X 3 )e l , where x t are the coordi¬ 
nates with respect to an orthonormal coordinate frame 
{e*} on D and {e 1 } denotes the standard orthonormal 
basis of Euclidean space. The spontaneous deformation 
tensor is given by 

Fij = (A - A - ") mrij + \~ v 5ij. (13) 

The reference metric determined by the nematic director 
field is 

9ij = {F )ik{F)kj 

= (A - A ~ u ) 2 n k n k ninj 
+ 2A '"'(A — A l ')nii r ij + A 2 u 8ki5kj- (14) 

Since rifcrifc = 1, we obtain 

9ij = (A 2 — A 2u )riinj + A 2v S^j. (15) 

Here we are interested in surface director field patterns 
on initially flat thin sheets, so we consider a plate of 
thickness h with a standard domain of parametrisation 
of the form H = wx ( — |, |) and director Helds of the 
form 

n{x) = ni(x\,X 2 )e 1 + n 2 (x i,x 2 )e 2 , (16) 

where w denotes the mid-plate. That is, we assume that 
the same surface director field is repeated at each level 
across the plate thickness. The reference metric now has 
components 

g a /3 = (A 2 - \~ 2v ) n a np + A ~ 2 v 8 a p, (17) 

and ^33 = A -2 ", g a 3 = 0. The reference fundamental 
forms in the reduced two-dimensional model are 

a a p{xi,x 2 ) = g a p(x-i, x 2 , 0 ), b a/ 3 (x 1 ,x 2 )= 0 . (18) 

Thus, in the thin sheet limit the components a a p of the 
first fundamental form of the mid-surface are given by 

«a /3 = (A 2 — \~ 2 u )n a np + X~ 2v S a p. (19) 


of a surface. In particular, it is uniquely determined by 
the first fundamental form a a p according to the equation 

k = Urlz - d 2 r 2 x + 

a 11 v 

- Tn r ? 2 + r? 2 r? 2 - r^r^). (21) 

A direct calculation of the Gaussian curvature associated 
with the nematic metric yields the expression 

K = i (A 2i/ - A -2 ) [(d$ip - d\ip - Adnpd 2 ip) sin(2^) 

+ 2 (did 2 '0 + (d 2 ip) 2 - (9i^) 2 ) cos(2ip)] , (22) 

in terms of the angle scalar field ip. 

Define the orthogonal dual to a director field to be 
the director field obtained by rotating all the nematic 
directors by 7r/2. That is, the dual director field is char¬ 
acterised by the angle field ip = ip + 7r/2. Upon tak¬ 
ing the orthogonal dual ip —> ip + 7t/ 2, we find that 
d a ip —» d a ip, sin(2'!/>) —» sin(2 -ip + tv) = — sin 2 ip and 
cos(2 ip) —)• cos(2 ip + 7 r) = — cos 2ip, so that 

K -+-K as ip ^ ip + tt/2. (23) 

That is, the orthogonal dual of any director field has 
precisely the opposite Gaussian curvature at every point. 

III. SHIFTED DIRECTOR PATTERNS 

Suppose that the components of the director field de¬ 
pend only on one of the coordinates, so that n = n( x 2 ), 
say. Such a director field can be specified as 

n = cosip(x 2 ) e 1 + sin^(a; 2 ) e 2 . (24) 

The Gaussian curvature of the associated metric is given 
by 

K = -l( A“ 2 - A 2 ") (ip" sin 2 ip + 2t// 2 cos 2ip) . (25) 

For constant I\ € R, this yields a second order ordinary 
differential equation in ip: 

d 2 

cos 2^ = 4C(K), (26) 


Since the 2D in-plane director field is modelled as a nor¬ 
malised vector field, it can be specified by a single angle 
scalar field ip = ip(x\,x 2 ) such that 


n\ = cosip(xi, x 2 ), n 2 = smip(xi,x 2 ). (20) 

Recall that the Gaussian curvature K at a point p on 
a surface is defined as the product of the two principal 
curvatures at p. In terms of the fundamental forms of the 
surface, it can be expressed as K = det (b a p)/ det(a aj g). 
According to the Theorem Egregium of Gauss, the Gaus¬ 
sian curvature is a characteristic of the intrinsic geometry 


where C(K) = K/(X~ 2 - A 21 '). This ODE is solved by 

ip(x 2 ) = ±^cos -1 (ci + c 2 x 2 + 2 C(K) xf ^, (27) 

where c\, c 2 are constants of integration. The corre¬ 
sponding director field generates constant Gaussian cur¬ 
vature K wherever it is well-defined. Fig. |T] shows ex¬ 
amples of director fields that generate constant Gaussian 
curvature of various types. 


Director fields of the form (24) can be generated by 


shifting a fixed curve along the aq-axis. We now develop a 



4 



FIG. 1. Director fields of the form n = cosip(x 2 ) e 1 + sin^(* 2 ) e 2 on the unit square for different choices of the function 
ip = tp(x 2 ). The resulting Gaussian curvature K upon stimulation is also indicated. 


more geometric description for director fields that can be 
generated by uniform translation of a curve 7 = 7 (f) = 
7 i(t)e 1 + 7 2 (f)e 2 along a fixed direction specified by a 
unit vector u = uie 1 + u 2 e 2 . For such a director field, 
it is natural to change coordinates from the Cartesian 
(xi, x 2 ) to (f, r ) where t is the parameter along the curve 
7 and r is the parameter in the direction of translation. 
That is, we have 


Xi(t,r) := 71 (f) + ru\ , x 2 (t, r) := 72(f) + ru 2 . (28) 


The director field n = nie 1 + n 2 e 2 at each point (f, r) is 
given by 


ni(t) = 


7i 




■il 


n 2 (t) = 


72 


VWT 


7 ? 


(29) 


Note that the independence of the director field com¬ 
ponents from r follows by construction. The metric 
components A = [a ai a] = [(A 2 — A ~ 2 v )n a np + A _ 2 "J Q/ 3 ] 
in Cartesian coordinates transform according to A —» 
J 1 AJ, where J is the Jacobian matrix 


3 = ( d * Xl J * 1 ) . (30) 

V d t x 2 o r x 2 J y J 

That is, 

( a tt a tr \ _ f 7i (t) u i \ T ( an ai2 \ ( 7i (t) ui \ 
\art a rr J \Y 2 (t) u 2 J \a 2 i a 22 J ^ 7 2 (t) u 2 J 

(31) 


Now for a given curve 7 and specified direction u , we 
can compute the Gaussian curvature as a scalar field K = 


K(t) using the equation 


K = — 


— (d t T r tr -d r r r tt ^ 

a t t v - 

_ r* r r 

1 tt L tr 


■ r‘ r r 

i- + 


T^r pr j^r 

1 tr L tr 1 tt L rr 


). (32) 


Here we note a remarkable result with an elegant ge¬ 
ometric interpretation. Take a circular arc 7 (t) = 

R (cost, sin t), where R = ^ 7 = (A -2 — A 2l/ ) 1,/2 , and K > 
0 is a constant. Now let u = (cos a, sin a), where a is 
some fixed angle. A direct computation of the Gaussian 
curvature associated with the director field generated by 
translating 7 along u , yields K exactly. It is assumed 
that the circular arc and direction of translation are cho¬ 
sen such that the director field does not self-intersect. 

Recall that the tractrix in the xix 2 -plane whose axis 
coincides with the 07 -axis is the planar curve 7 passing 
through a point ( 0 , a) with the property that the length 
of the segment of the tangent line from any point on the 
curve to the uq-axis is constant and equal to a > 0. One 
parametric representation of the tractrix is provided by 


71 (f) = a (tanht — t), 72 (f) = asechf. (33) 


We now consider the director pattern that is generated 
by translating a tractrix along its axis. That is, we 
take 7 to be as in equation (33) and choose the di¬ 
rection of translation to be « = ( 1 , 0 ) as shown in 
Fig- U A direct computation of the Gaussian curvature 
of such a pattern generated by a tractrix with parameter 
a = (^ -2 — A 2y ) 1//2 with K < 0, yields I\ exactly. 
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a simultaneous expansion in the orthogonal directions. 
Thus, one might expect that the intrinsic curvature prop¬ 
erties of a patterned sheet will be qualitatively reversed 
when the pattern is replaced by its orthogonal dual, even 
though the scaling factors A and A -v are not exact re¬ 
ciprocals when v ^ 1. Indeed, we have shown that this 
reversal is exact in the sense that the Gaussian curva¬ 
ture distribution becomes precisely the negative of the 
dual pattern. 

It is also interesting to connect these ideas to the me¬ 
chanical response of a patterned nematic sheet that is 
being cooled instead of heated. On cooling an initially 
flat nematic sheet, the rod-like molecules expand along 
the director and contract in the orthogonal directions. So 
once again one may expect the activated surface to ex¬ 
hibit the opposite of the curvature properties that mani¬ 
fest on heating. Quantifying this opposite response how¬ 
ever is somewhat delicate, since the scaling coefficients 
are typically temperature dependent and cooling a flat 
sheet by a particular temperature will not necessarily 
produce the exact opposite of the response observed by 
heating the flat sheet by the same temperature. 


FIG. 2. The nematic pattern obtained by shifting a cir¬ 
cular arc along a fixed direction generates constant positive 
Gaussian curvature K > 0. The nematic pattern obtained by 
shifting a tractrix along its axis generates constant negative 
Gaussian curvature K < 0. 


IV. ORTHOGONAL DUALITY 


The result that the director fields in Fig. [2] generate op¬ 
posite Gaussian curvature may seem surprising on first 
inspection, since the patterns look somewhat similar. 
However, in light of the orthogonal duality result (231, 


we find that the patterns make perfect sense. In partic¬ 
ular, the positive curvature pattern is equivalent to the 
pattern obtained by rotating the sheet by 90°. One can 
see that this equivalent pattern is exactly the orthogo¬ 
nal dual of the negative curvature pattern of Fig. [2j as 
illustrated in Fig [3] 



FIG. 3. Orthogonal duality of director fields. The direc¬ 
tor fields that were discovered to generate constant Gaussian 
curvatures of opposite sign are seen to be equivalent to or¬ 
thogonally dual patterns. Note that every tractrix intersects 
every circle at right angles. 


This observation is consistent with physical intuition 
since the mechanical response of nematic solids to stim¬ 
uli such as heat is a contraction along the director and 


V. WILLMORE FUNCTIONAL 

It is well known that for a surface in R 3 , the compo¬ 
nents a a p and b a p of the first and second fundamental 
forms satisfy the Gauss-Codazzi-Mainardi equations 

dpT aaT — d a r a p T + 

— botcrbp T b a pb(j T , 

dpb aa - d a b a p + T^bpfj. - T = 0, (34) 

where V a p T = §( dpa aT + d a ap T - d T a a p) and T a ap = 
a aT T a p T . Furthermore, any pair (a, b) consisting of a 
symmetric and positive definite matrix field ( a a p ) and 
a symmetric matrix field ( b a p ) that satisfy the Gauss- 
Codazzi-Mainardi equations determines a unique surface 
up to a rigid transformation in R 3 mi. 

To determine the equilibrium configuration of the mid¬ 
surface of an initially flat nematic sheet upon stimulation, 
we need to know the components b a p of the second fun¬ 
damental form that minimise the bending energy 

£[b] = \ [ A^bcpb^dS (35) 

d Ju 

subject to the Gauss-Codazzi-Mainardi constraints de¬ 
fined by the metric (a a p) generated by the director field. 
This functional can be rewritten in terms of the mean 
and Gaussian curvatures of the mid-surface as 

e =lL{ ^- 2K ) dS ‘ <36) 

Since a a p = a a p, and the Gaussian curvature K is an 
isometric invariant, the problem reduces to minimising 
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the Willmore functional 


written as 


iw 


[ H 2 dS 

(37) 

|-27t rR sinh ^^/\K\pj 

J UJ 


Iw ~ Jo Jo vm 


among isometric immersions of the metric 

The Willmore functional can be written in terms of the 
principal curvatures K \, k 2 as 

Iw — j (k i + 1^2) 2 dS 

^ J UJ 

f {Ki~K 2 ) 2 dS+ ( KdS, (38) 

J UJ J UJ 


1 

4 


(k 2 + « 2 ) dpd6. (40) 


It is well-known that at every point of a hyperbolic 
surface, there exists a pair of asymptotic curves along 
which the normal curvature vanishes m- Furthermore, 
the isometric immersions of the hyperbolic metric of con¬ 
stant Gaussian curvature K are in correspondence with 
solutions of the sine-Gordon equation: 


d 2 f> 

dudv 


= \K\ siiuf>(u, v), 


(41) 


where K = KiK 2 is the Gaussian curvature. We note 
that if the condition k\ = k 2 throughout to is consistent 
with the metric, then the equilibrium configuration 
corresponding to minimal bending energy is determined 
by imposing this condition. The sphere is the only 
compact surface in R 3 whose principal curvatures are 
equal everywhere [16]. Since surfaces of the same 
constant Gaussian curvature K are locally isometric 
by Minding’s theorem a flat nematic sheet whose 
director field encodes constant positive curvature K 
will form part of a sphere of radius R = 1 /\[K upon 
stimulation, assuming that the sheet is small enough to 
exclude the possibility of self-intersection as in Fig. [4] 


/ 



FIG. 4. An initially flat nematic sheet forming part of a sphere 
upon stimulation. 

In the case of a thin sheet encoded with constant nega¬ 
tive curvature K < 0, the study of the equilibrium shapes 
demands a considerably more delicate analysis. We fol¬ 
low Gemmer et al. 03 and consider equilibrium con¬ 
figurations of hyperbolic elastic disks of radius R that 
have already undergone local deformations. The equilib¬ 
rium configurations in the thin sheet limit correspond to 
minimisers of the Willmore functional or, equivalently, to 
minimisers of 


iw — / ( K i + K !) ds, 

J UJ 


(39) 


over isometric immersions of the metric. In geodesic po¬ 
lar coordinates (p,9), where the metric takes the form 

ds 2 = dp 2 + \K\ 1 sinh 2 (^J\K\p s j d6 2 , Eq. (391 can be 


where u, v are coordinates for a local parametrisation by 
asymptotic curves and (f> denotes the angle between the 
asymptotic curves. The angle <f> is related to the principal 
curvatures by 

«i = \K\ tan 2 r\ = \K\ cot 2 (42) 


so that 


r 2 -K r R sinh 


Iw = 


(vw\p) 


VW\ 


tan — + cot — j dpdO. 

(43) 


In [18] , the minimisation of Eq. (43) over all smooth 
solutions to the sine-Gordon equation is analysed numer¬ 
ically to conclude that the principal curvatures of smooth 
isometric immersions of a hyperbolic disk of radius R sat¬ 
isfy 


max(M,M)> (44) 

64 

Moreover, it is shown that this lower bound is attained 
by surfaces that are geodesic disks lying on hyperboloids 
of revolution of constant Gaussian curvature K. 

Hyperbolic surfaces whose two asymptotic curves are 
straight lines that intersect at an angle 9 form a one- 
parameter family of surfaces Ag known as Amsler sur¬ 
faces [20'. Each such surface is uniquely determined by 
the angle 9 between the asymptotic lines. In [18] , the 
odd periodic extension of a subset of the Amsler surface 
An/n bounded between the asymptotic lines is taken to 
generate a periodic profile with n waves. The resulting 
n-periodic shapes are referred to as periodic Amsler sur¬ 
faces A ni which are not smooth as they have discontinu¬ 
ities in their second derivatives at the lines of inflection. 

Numerical investigations of the bending energy associ¬ 
ated with geodesic disc cut outs of A n have revealed that 
for large values of e = \J\K\R 1 particular ?r-wave pe¬ 
riodic Amsler surfaces are energetically more favorable 
than the smooth saddle shapes that correspond to discs 
lying on hyperboloids of revolution. Furthermore, it has 
been shown that there exists a sequence of critical radii 
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R n ~ ln(n) beyond which the principal curvatures and 
hence Willmore functional of A n diverge. If the radius of 
the hyperbolic disc is sufficiently large for periodic Am- 
sler surfaces to be energetically more favorable than the 
hyperboloid solutions, the particular n -wave surface that 
is expected to form is predicted by identifying the small¬ 
est critical radius R n greater than R. Thus, as the ra¬ 
dius of the hyperbolic disc increases one observes that 
n-periodic shapes with an increasing number of waves 
become more energetically favorable mm- 

It is interesting to contrast the case of shape selection 
of sheets encoded with constant positive curvature with 
that of sheets patterned with a hyperbolic metric. In 
the constant positive curvature case, the analysis sug¬ 
gests that subsets of spheres are the most energetically 
favorable isometric immersions provided that the dimen¬ 
sions of the initial sheet are small enough to rule out 
self-intersection of the deformed surface. That is, the 
nature of the solution is not influenced by the shape and 
size of the initially flat sheet. In the hyperbolic case, 
however, we notice that the extrinsic geometry of the de¬ 
formed surface is crucially dependent on the particular 
shape and size of the initially flat sheet. In particular, 
surfaces with dramatically different extrinsic geometries 
may form by simply increasing the size of the sample 
patterned with a particular intrinsic geometry. 


VI. BLUEPRINTING NEMATIC SHEETS OF 
PRESCRIBED GAUSSIAN CURVATURE 

One can use the shifted circular arc patterns of Fig. [2] 
to blueprint a metric of constant positive Gaussian 
curvature that is realised upon stimulation of a thin 
nematic sheet. Since the magnitude of the Gaussian 
curvature is determined by the radius of the shifted 
circular arcs, we need to use a shifted pattern consisting 
of circular arcs that are C-smoothly joined in order to 
extend the pattern to cover larger domains, as shown in 

Fig. |5] 

The pattern on the larger domain encodes positive 
Gaussian curvature at all points except along the lines 
where the distinct circular arcs join. At these points, the 
Gaussian curvature is not defined and sharp creases may 
develop upon stimulation. Similarly, one can extend the 
shifted tractrix patterns of Fig. [2]to encode constant neg¬ 
ative curvature on extended domains as shown in Fig. [6] 
As in the positive curvature case, the Gaussian curvature 
is not well-defined along the solid horizontal lines indi¬ 
cated in Fig. [6] as the metric components are not contin¬ 
uously differentiable at these points. 

In the case of the circular arc patterns, for instance, 
note that the unit tangents to the arcs give the director 
field whose first derivative at each point coincides with 
the vector pointing towards the centre of the circular arc 
at that point and so clearly exhibits a jump discontinuity 
across the solid horizontal line. The jump discontinuities 


O 



FIG. 5. Blueprinting prescribed positive Gaussian curvature 
K > 0 on a rectangular sheet of arbitrary dimensions whose 
edges are indicated by the dashed lines. The Gaussian cur¬ 
vature is well-defined everywhere except along the solid hori¬ 
zontal lines where the distinct arcs join. 


of the extended hyperbolic pattern follow by orthogonal 
duality. 



FIG. 6. Blueprinting prescribed negative Gaussian curvature 
K < 0 on a rectangular sheet of arbitrary dimensions whose 
edges are indicated by the dashed lines. The Gaussian cur¬ 
vature is well-defined everywhere except along the solid hori¬ 
zontal lines. 


VII. NEMATIC PATTERNS ON CURVED 
SURFACES 

Often we are interested in shape transitions of sur¬ 
faces that are initially curved. In such cases the surface 
may have a non-flat first fundamental form at the out¬ 
set, which would then change due to the activation of 
some local deformation pattern. In the context of thin 
nematic shells, we would start with a shell whose mid¬ 
surface is given as an embedding in three dimensions. 
This embedding determines the first fundamental form 
as a non-trivial (0, 2) metric tensor, which would then 
change to a new reference metric a a p that depends on 
the director field pattern on the surface. 

Here we highlight an interesting result on the effect of 




















two sets of director fields on the Gaussian curvature of 
surfaces of revolution. Consider a surface of revolution 
whose standard parametrisation is given by 


r(u,v) = {f(v)cosu,f(v)sinu,g(v)), (45) 


where the pair of real-valued functions (/, g) specifies the 
profile curve of the surface. The associated metric is 
given by 

ds 2 = r u ■ r u du 2 + 2 r u ■ r v dudv + r v ■ r v dv 2 

= f 2 du 2 + (f' 2 +g' 2 )dv 2 . (46) 

A direct calculation of the Gaussian curvature yields 


f'g'g" - f"g ' 2 

f{f 2 +g ' 2 ) 2 ' 


(47) 


Now suppose that the surface is patterned such that the 
directors are aligned with the v = constant lines. That 
is, the director field lines are circles centred on the axis 
of revolution. The metric upon stimulation becomes 


ds 2 = f 2 X 2 du 2 + (, f 2 + g' 2 ) X~ 2u dv 2 , 


and the Gaussian curvature transforms as 


f'g'g” - f"g ' 2 
' f(f ' 2 + g' 2 ) 2 ' 


(48) 


(49) 


That is, for any surface of revolution patterned by di¬ 
rector fields aligned with the horizontal circles centred 
on the axis of revolution, the Gaussian curvature trans¬ 
forms according to the elegant formula K —► A 2 "A'. This 
result was established in [223 for the special case of a 
sphere patterned by azimuthal director fields. Here we 
see that the result is valid for all surfaces of revolution, 
including surfaces of negative and non-constant Gaussian 
curvature. 

For the orthogonal dual of the considered pattern on 
any given surface of revolution, the metric upon stimula¬ 
tion is given by 

ds 2 = f 2 X~ 2v du 2 + (/' 2 + g' 2 ) A 2 dv 2 , (50) 


and a direct calculation of the Gaussian curvature yields 


K = A -2 


f'g'g" - f”g ' 2 

f(f ' 2 + g ' 2 ) 2 ‘ 


(51) 


That is, for any surface of revolution patterned by direc¬ 
tor fields aligned with the profile curves of the surface 
(i.e. lines of constant u), the Gaussian curvature trans¬ 
forms according to K —» X~ 2 K. Note the remarkable 
fact that this transformation is independent of the opto- 
thermal Poisson ratio v. In the case of a sphere, such 
patterns correspond to the director fields being aligned 
with the lines of longitude. 

In the special case of a sphere of radius R , these results 
show that upon the stimulation of an azimuthal pattern 
the Gaussian curvature decreases to A 2 "/!? 2 , whereas for 


the orthogonal dual pattern the Gaussian curvature in¬ 
creases to A ~ 2 /R 2 . We now show how these two patterns 
and the associated Gaussian curvature responses can be 
viewed as the extreme cases of the family of loxodromic 
spiral patterns on a sphere, as shown in Fig. [7] 

Recall that a loxodrome on a sphere is the curve that 
cuts a meridian at a constant angle a. If a is not a right 
angle, then the loxodromes form spirals on the sphere. 
We assume that a sphere of radius R is patterned with a 
director field whose integral curves consist of loxodromes 
of angle a , so that n = cos a eg + sin a e^, where 9 and 
cf> are the polar and azimuthal angles, respectively, and 
eg, e ^ are the corresponding unit tangent vectors in R 3 . 
Upon stimulation, the point ( 9 , <j >) transforms according 
to (9, <p) T —> F(9 , (f>) T , where F is the matrix representa¬ 
tion of the spontaneous deformation tensor with respect 
to the orthonormal basis eg, e^. The round spherical 
metric ds 2 = R 2 (d9 2 + sin 2 9 d<f> 2 ) transforms as 

ds 2 ( d 0 d 0 ) F T aoF ^ = (d 0 dcf) a (^j , (52) 
where 


0 

R 2 sin 2 9 


(53) 


Directly computing the Gaussian curvature using the new 
metric a = F T a$F, we obtain the formula 




2 


cos 2a 


1 

R 2 ' 


(54) 


We notice that the resulting Gaussian curvature is 
a constant between the values A 2v /R 2 and A ~ 2 /R 2 
depending on the angle a. Note that although the 
Gaussian curvature upon stimulation is constant and 
positive in each case, the resulting surface is generally 
not a sphere. For example, in the case of azimuthal 
patterns (a = 7r/2) a spindle or “thorny sphere” of the 
corresponding curvature is expected to form : 22 j. 


VIII. COMMENTS 

It is often the case that key physical characteristics 
of materials, such as adhesive and optical properties, 
are determined by surface structure. The engineering of 
switchable surfaces consisting of active materials such as 
nematic glasses offers the possibility of controllably and 
reversibly changing the surface geometry of thin struc¬ 
tures for use in a variety of potential devices. Possible 
applications may include microfluidic mixers and pumps, 
adjustable optical lenses and switchable textured surfaces 
for use in tablet computers. There are two key mecha¬ 
nisms for generating shape transformations in thin sheets 
of patterned nematic glasses. One method is by means 
of defects in the director field pattern and the other is 
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FIG. 7. Loxodromic spiral director fields of the form n = cos a eg+sina on a sphere of radius R for the cases a = 0, 7r/4, 7t/2 
and the Gaussian curvature generated upon stimulation. Note that a is the angle that the directors make with the meridian 
lines. 


by using a smooth in-plane director field that generates 
Gaussian curvature. 

In this paper, we have clearly highlighted smooth di¬ 
rector fields that generate constant Gaussian curvature 
of any desired value, with the hope that such patterns of¬ 
fer feasible initial targets for experimentalists to confirm 
the predictions of the theory and test its limitations. We 
have also identified an important symmetry in the Gaus¬ 
sian curvature associated with a given nematic pattern 
on an initially flat sheet, namely that the orthogonal dual 
of a given pattern generates the exact opposite Gaussian 
curvature at every point. This observation has an impor¬ 
tant practical implication for experimentalists, who can 
produce a sheet whose Gaussian curvature distribution 
is the exact opposite of a given patterned sheet by sim¬ 
ply adding n/2 radians to the nematic alignment angle 
associated with the given sheet. If the viability of pro¬ 


ducing switchable nematic surfaces using such patterns 
is confirmed, one may begin to effectively combine more 
elaborate smooth director fields and topological defects 
to devise a variety of powerful shaping mechanisms for 
use in a wide range of potential applications. 
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